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Abstract. In this paper, we study the feasibility of conducting quantum key 
distribution (QKD) together with classical communication through the same optical 
fiber by employing dense-wavelength-division-multiplexing (DWDM) technology at 
telecom wavelength. The impact of the classical channels to the quantum channel 
has been investigated for both QKD based on single photon detection and QKD based 
on homodyne detection. Our studies show that the latter can tolerate a much higher 
level of contamination from the classical channels than the former. This is because 
the local oscillator used in the homodyne detector acts as a "mode selector" which can 
suppress noise photons effectively. We have performed simulations based on both the 
decoy BB84 QKD protocol and the Gaussian modulated coherent state (GMCS) QKD 
protocol. While the former cannot tolerate even one classical channel (with a power 
of OdBm), the latter can be multiplexed with 38 classical channels (QdBm power per 
channel) and still has a secure distance around lOfcm. Preliminary experiment has 
been conducted based on a lOOA4Hz bandwidth homodyne detector. 
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1. Introduction 

One important practical application of quantum information is quantum key distribution 
(QKD), whose unconditional security is based on the fundamental laws of quantum 
mechanics [UEj- Today most of the fiber-based QKD experiments are conducted through 
dedicated "dark" fibers. As the applications of QKD have been extended from point 
to point links to network configurations [3j HJ |5], it would be much more appealing 
if QKD can be conducted through the existing fiber network together with classical 
communication signals. In this "coexistence" architecture, the quantum signal for QKD 
shares a common optical fiber with unrelated classical signals [6]. 

The coexistence architecture based on wavelength-division-multiplexing (WDM) 
technology was proposed by Townsend in the late nineties [7] and has been studied 
thereafter [HI El EJ QUI ITT] . These studies show that the existence of strong classical 
traffic could be detrimental to the QKD channel. This is because a classical signal is 
typically many orders of magnitude stronger than the quantum signal. Thus even a 
small crosstalk from a classical channel could overwhelm normal QKD operation. One 
feasible way is to place the quantum signal at the 1300nm "original" or O-band while 
placing the classical signal at the 1550nm "conventional" or C-band [6j. For such a 
large wavelength separation, both the leakage of the classical signal and the anti-stokes 
Raman scattering can be effectively suppressed to a tolerable level. 

However, in many cases there are advantages to place both the quantum and the 
classical signals in C-band. For example, the fiber loss at C-band is significantly lower 
than that at O-band, so the secure distance of QKD could be extended. Furthermore, 
this architecture is more compatible with today's fiber network. 

In [in], a BB84 QKD system is multiplexed with classical signals using C- 
band reconfigurable optical add drop multiplexer (ROADM). In [TTJ, the authors 
multiplexed four classical channels with one quantum channel using C-band 100GHz 
dense- wavelength-division- multiplexing (DWDM). This is a special case where the four 
classical channels are used by the legitimate users (Alice and Bob) for key distillation and 
encrypted communication. The total power of the four classical channels is only about 
—22dBm. By placing the quantum signal at a non-adjacent channel of the classical 
signal, the maximum secure distance has been shown to be around 40km [UJ. These 
previous studies suggest that with standard technology, it is still infeasible to place both 
QKD signal and strong classical signal (OdBm) at C-band. 

We remark that all the previous studies have been conducted in QKD systems 
employing single photon detector (SPD), implementing, for instance, the standard BB84 
QKD protocol. Here, we will show that QKD systems based on optical homodyne 
detection are inherently robust against noise due to multiplexing. Intuitively, this is 
because a strong local oscillator (LO) is used to beat with the weak quantum signal 
during homodyne detection. Only photons in the same spatiotemporal and polarization 
mode as the LO can be detected while noise photons in different modes will be suppressed 
effectively. Thus the LO acts as a "mode selector" [12]. We remark that in a recent 



Feasibility of QKD through DWDM network 



3 



free space QKD experiment, an optical homodyne detector was employed to suppress 
ambient light [T3] . 

In this paper, we study quantitatively the amount of noise photon added into the 
QKD channel in a typical DWDM setup. All the theoretical simulations are carried 
out based on the performance of commercial components. Two specific QKD protocols, 
namely, the decoy state BB84 protocol [HI EES! EE] and the Gaussian modulated coherent 
state (GMCS) QKD protocol [17] have been chosen for numerical simulations. The BB84 
QKD protocol [1] is the first and the most well known QKD protocol. Its unconditional 
security has been rigorously proved based on the laws of quantum mechanics [18], even 
when implemented on practical setups with some imperfections, such as weak coherent 
pulse, detector dark counts [191 I2Q] and detector efficiency mismatch [2TJ . The security 
of the GMCS QKD protocol was first proven against individual attacks with direct 
|22] or reverse (TTJ, |23] reconciliation schemes. Security proofs were then given against 
general collective attacks [2U [25j [26] . To date, three groups have independently claimed 
that they have proved the unconditional security of GMCS QKD [27j . 

Our simulation results show that it is possible to multiplex the GMCS QKD 
with a OdBm classical channel using a C-band 100GHz DWDM without significantly 
reducing its performance. Even multiplexed with 38 classical channels (OdBm power 
each channel), the GMCS QKD still has a secure distance around 10km. Preliminary 
experiment has been conducted based on a 100MHz bandwidth homodyne detector. 

We remark that although our studies are conducted in the GMCS QKD, our 
conclusions should be applicable in other QKD protocols employing homodyne detection 
such as the discrete modulated continuous variable QKD protocol [28] . 

This paper is organized as follows: Section 2 contains theoretical analysis of various 
noise sources in a generic quantum/classcial signals coexistence scheme. In Section 3, 
We will quantify the contribution from each noise source based on the performance 
of commercial components. In Section 4, we will present simulation results in both 
the decoy state BB84 protocol and the GMCS QKD protocol. We will also present 
preliminary experimental results with a 100MHz shot noise limited homodyne detector. 
Section 5 is a brief conclusion. 



2. Theoretical analysis on noise contribution 

In a typical coexistence architecture based on DWDM, the noise photons in quantum 
channel can be contributed by several sources [HI EEl EH] , including the leakage photons 
from the classical channels due to the finite isolation of the DWDM components, the 
"in-band" noise photons generated in optical fiber from nonlinear processes, such as 
four-wave mixing (FWM) and the spontaneous Raman scattering, the in-band amplified 
spontaneous emission (ASE) photons generated by optical amplifiers. Here, in-band 
noise refers to noise photons within the spectral bandwidth allocated to the quantum 

| A typical 40-channel 100GHz DWDM system with 1 channel for quantum communication and 38 
nonadjacent channels for classical communication. 
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Figure 1. A typical scheme for multiplexing quantum channels with classical channels. 
MUX- multiplexer; DEMUX-dcmultiplexer; EDFA-erbium-doped fiber amplifier. 



channel. In this section, we will quantify the amount of noise photons contributed 
by each of the above sources based on a typical DWDM configuration as shown in 
Fig.l. In Fig.l, we assume that an erbium-doped fiber amplifier (EDFA) is employed to 
boost the optical power of classical channels before multiplexing with quantum channels. 
Furthermore, we assume that all the classical channels are placed at wavelengths longer 
than that of the quantum channels, since the spontaneous anti-Stokes Raman scattering 
(SASRS) is typically weaker than the spontaneous Stokes Raman scattering. 

In this paper, we assume that the eavesdropper (Eve) can control all the classical 
channels and the EDFA (see Fig.l) but she cannot access the multiplexer (MUX) and 
the demultiplexer (DEMUX) used for multiplexing the quantum signals with classical 
signals. One special example is that the classical signals are actually used by Alice 
and Bob for authentication, error correction and privacy amplification [IT]. In the 
more general cases where the classical channels are allocated to other users, Alice and 
Bob can place the MUX and DEMUX in their local secure stations. Under the above 
assumptions, the quantum signals sent by Alice are calibrated after the MUX (point A 
in Fig.l). So the performance of the QKD system is independent of the insertion loss 
of MUX. On Bob's side, the insertion loss of the DEMUX can be treated as part of the 
loss in Bob's detection system. 

2.1. Amplified spontaneous emissions (ASE) of EDFA 

It is well known that an ideal, noise-free amplifier cannot exist [29J. In the case of an 
optical amplifier, the fundamental noise originates from the spontaneous emission. The 
ASE from a practical EDFA has a broad bandwidth on the order of tens of nm which 
can be treated as a broadband noise source with a flat spectral power density within the 
spectral bandwidth of the quantum channel. We remark that a practical laser source 
also has a broadband noise background, which can be modeled as ASE from a virtual 
optical amplifier. 

The average ASE photon number in one spatiotemporal mode is given by [30] 

(N ASE ) = 2n sp (G - 1) (1) 

Here the factor 2 accounts for the two orthogonal polarization modes. G is the gain of 
the EDFA, n sp > 1 is the spontaneous emission factor. If the spontaneous emission is 
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the only noise source (no excess noise), n sp = 1. 

In practice, the excess noise of an EDFA is commonly quantified by its noise figure 
(NF). In the unsaturated regime, NF is related to n sp by [30] 

NF = 1 + 2 "- (G - 1) (2) 
G 

In the high gain range (G 3> 1), NF ~ 2n sp . 

Typically, the ASE power is much lower than that of the classical signal. However, 
its bandwidth is much broader and extends into the quantum channel. Thus the ASE 
will contribute to in-band noise. Fortunately, the MUX used at Alice's side functions 
as a bandpass filter and can greatly suppress this in-band ASE noise. Given the cross 
channel isolation of the MUX is £i, the in-band ASE photon number (per spatiotemporal 
mode) measured at the output of the MUX (point A in Fig.l) is 

(N<$ E ) = 2&n v {G - 1) (3) 

Note, in this paper, we will not consider the "out-of-band" ASE noise photons, 
since they are typically much weaker than the classical signals themselves. 



2.2. Leakage from classical channel 

Although the classical signal has a different wavelength as the quantum signal, a small 
fraction of the classical signal will leak into the quantum channel due to the finite 
isolation of the DEMUX. In a BB84 QKD system, this leakage will contribute to out- 
of-band noise, which could be further reduced by using spectral filters at the receiver's 
end. In a GMCS QKD system, this leakage contributes noise photons in "unmatched 
mode" of the LO. 

We define the power of the classical signal output from the communication fiber as 
P out (measured at point B in Fig.l). Given the isolation of the DEMUX £ 2 , the power 
of the leakage signal received by Bob (measured at point C in Fig.l) is Pi ea k = £2 -Pout- 
The average leakage photon number per second is 

<iCi> = (4) 
Here h is Planck's constant, and v is the frequency of the classical signal. 



2.3. Spontaneous anti-Stokes Raman scattering (SASRS) 

As the strong classical signals propagate along the optical fiber, noise photons at 
different wavelength can be generated through various nonlinear optical processes. If 
the wavelength of the noise photons coincides with that of the quantum signal, they 
cannot be filtered out at the receiver's end and will contribute to in-band noise. It has 
been shown that SASRS is the dominant nonlinear process when the quantum channel 
is placed at the shorter wavelength of the classical channel [6j [11] . 
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The SASRS noise power within a bandwidth of A A (measured at point B in Fig.l) 
is given by [6] 

Psasrs = Pinf3zr] ch A\ = P out f3zAX (5) 

Here (3 is the spontaneous Raman scattering coefficient, Pj„ (measured at point A in 
Fig.l) is the input power of the classical signal, z is the fiber length and r] ch is the 
transmittance of the optical fiber. 

To estimate the noise photon number per spatiotemporal mode, we first use the 
relation v — c/X to determine the total mode number corresponding to a bandwidth of 
AA and a time window of At = Is to be N moc i e = \AuAt\ = p-AA. Here c is the speed 
of light in vacuum. 

Given the insertion loss of the DEMUX is tjdmu, the in-band SASRS photon number 
(per spatiotemporal mode) measured at the output of the DEMUX (point C in Fig.l) 
can be calculated from (5): 

/7\r(c) \ Psasrs A 3 

( n sasrs) = huN Vdmu = PoutP ZTlDMU (6) 
Again, in the derivation of (6), we have used the relation v = c/X. 

2.4- Four-wave mixing 

Four-wave mixing (FWM) is a third order nonlinear process generated by the x 
nonlinearity of the optical fiber when two or more pumps exist. For FWM process 
to be efficient, phase-matching condition is required. Although FWM could be the 
major noise source at very short distance, it is much weaker than Raman Scattering 
for a practical fiber length [JO]. Furthermore, FWM can be effectively suppressed by 
optimizing the channel configuration [TO] ITT] or using polarization multiplexing [10]. In 
this paper, we simply neglect FWM. 



3. Experimental characterization of various noise sources 

3.1. The performance of commercial MUX/DEMUX 

Theoretical studies in Section 2 show that the noise level is dependent on the 
performance of MUX/DEMUX: cross-channel isolation and insertion loss. 

We have tested two commercial C-band 100GHz MUX/DEMUX from JDSU 
(model number: WD1508D1B). Each of them has 8 channels with a channel separation of 
0.8nm (or 100GHz). The 3dB bandwidth of each channel is around 0.6nm (or 75GHz). 
Based on the availability, a 1559. 79nm fiber optic source module (ILX Lightwave, 
79800D) has been used as the laser source for the classical channel. Thus, we allocate 
channel 2 of the MUX/DEMUX as the classical channel and channel 8 as the quantum 
channel. The isolation £i of MUX is determined by sending a calibrated laser beam 
(with a wavelength of A 8 ) into channel 2 and measuring the optical power output from 
its common port. Similarly, the isolation £2 of the DEMUX is determined by sending 
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Figure 2. Experimental setups. Att-tunable optical attenuator; SMF-single mode 
fiber. 



a calibrated laser beam (with a wavelength of A2) into its common port and measuring 
the optical power output from its channel 8. In Table 1, we list the central wavelengths 
(A), insertion losses (L) and the cross-channel isolations (£) of two channels to be used 
in our experiments. 



3.2. Noise photons contributed by amplified spontaneous emissions of EDFA 

The noise contributed by an EDFA has been studied in Section 2.1. However, the level 
of noise photon given in (3) is too low to be measured with a conventional optical power 
meter directly. Instead, we have performed an experiment based on a modified setup to 
test the validity of (1). The modified experimental setup is shown in Fig. 2a. The EDFA 
is a commercial low noise fiber amplifier with a NF of 5.5dB (PriTel, LNHPFA-30-M). 
A tunable optical attenuator (Att. in Fig. 2) is used to simulate the transmission loss 
experienced by the classical signal before it reaches the multiplexer. 

Referring to Fig.2a, by setting the gain of EDFA to be G = 100, we have measured 
the noise power output from channel 8 to be P^se = —25dBm. On the other hand, 
from (1) the ASE photon is determined to be 351 per spatiotemporal mode. The energy 
of a 1550nm photon is hv = 1.28 x 10~ 19 J. So, the expected noise power within 75GHz 
bandwidth (corresponding to 0.6nm) is 351 x 75 x 10 9 x 1.28 x 10~ 19 J = —24.7dBm. 
Considering that the insertion loss of channel 8 is 0.9dB, the expected noise power 
Pase"^ = —25.6dBm, which reasonably matches with the experimental result P^e ■ 

Next, we determine whether the background noise of the laser source itself will make 
a significant contribution after being amplified by the EDFA. We first determined the 
background noise of the laser source. By using a tunable optical attenuator, the power 
of the classical signal (input to the EDFA) was set to —20dBm. The laser background 
noise at 1554. 94nm (As) within a 0.6nm wavelength range (corresponding to the channel 
bandwidth of MUX/DEMUX) has been determined to be —70dBm. The number of in- 
band laser noise photon in one spatiotemporal mode can be determined to be around 
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0.01. If the gain of EDFA is G, we expect that the number of the amplified laser noise 
photon is around 0.01G per spatiotemporal mode, which is significantly smaller than 
the number of ASE noise photon (see (1)). In this paper, we will simply neglect the 
contribution of the amplified laser noise. 

3.3. Spontaneous anti-Stokes Raman scattering generated in SMF28 fiber 

The SASRS noise photon number can be calculated from (6). The spontaneous 
Raman scattering coefficient of standard SMF28 fiber has been determined using 
the experimental setup shown in Fig.2b. The laser wavelength is 1559. 79nm and the 
output power after the MUX (point A in Fig. 2b) was set to AdBm. We measured the 
output power from channel 8 of the DEMUX for different fiber links: a 20km SMF28 
fiber spool and a 40km SMF28 fiber spool. Using (5), the spontaneous Raman scattering 
coefficient has been determined to be 2.85 x I0~ 9 (km ■ nm)^ 1 . Our result matches with 
the results reported in [UJ, which is between 2 x 10 -9 ~ 4 x 10~ 9 (km ■ nm)~ l depending 
on the Raman shift. 

In Section 4, we will calculate the secure key rates of both the decoy state BB84 
protocol and the GMCS QKD protocol in a typical DWDM configuration as shown in 
Fig.l. The simulation parameters are summarized in Table 2. To make our simulation 
results more applicable, some parameters in Table 2, such as NF, t]mux, Vdmu, £i 
and £2, have been chosen to be slightly worse than the experimental values on purpose. 
Since the spontaneous Raman scattering coefficient /3 is wavelength dependent, we have 
assumed the worst case of = 4 x I0~ 9 (km ■ nm)~ l [UJ. The parameters of the GMCS 
QKD system are from [32J. 

4. System comparison: single-photon detection scheme vs. homodyne 
detection scheme 

4-1. A single-photon detection based scheme: decoy state BB84 QKD with a weak 
coherent source 

Refer to Fig.l, in the BB84 QKD protocol, the quantum signal is at the single photon 
level, so the crosstalk between the quantum channels is negligible. We consider the 
simplest case where only one classical channel (at a longer wavelength and non-adjacent 
channel) is multiplexed with quantum channels. 

In the BB84 QKD, single photon detectors are employed to detect quantum signals. 
At telecom wavelength, InGaAs APDs working at the gated Geiger mode are frequently 
used as SPDs. In this case, the gating window of the SPD functions as a temporal filter, 
which can reduce the effective noise photon number. For other non-gated detection 
systems, this noise reduction can be achieved by introducing an adjustable detection 
time window. 

Given the channel bandwidth Au of MUX/DEMUX and the gating window At of 
the SPD, the total mode number of the noise photons which can be detected by the 



Feasibility of QKD through DWDM network 
Table 2. Simulation Parameters 



9 



Parameter 



Value 



NF (Noise figure of EDFA) 

a (Fiber attenuation coefficient) 

(3 (Spontaneous Raman scattering coefficient) 

Vmux (Transmittance of MUX) 

Vdmu (Transmittance of DEMUX) 

(Isolation of MUX) 
£ 2 (Isolation of DEMUX) 
Av (3dB channel bandwidth) 
At (Gating window of SPD) 
Va (Modulation variance, GMCS) 
VBob (Transmittance of Bob's system, GMCS) 
e (GMCS parameter) 
v e i (GMCS parameter) 

7 (Efficiency of reverse reconciliation algorithm) 



4 (or 6dB) 

0.21dB/km 

4 x 1CT 9 / (km ■ nm) 

0.71 (or 1.5dB loss) 

0.71 (or 1.5dB loss) 

10~ 8 (or -80dB) 

10~ 8 (or -80dB) 

75GHz 

Ins 

10 

0.6 

0.01 

0.01 

0.9 



SPD is 

N Mod = AuAt (7) 

The number of noise photons arrived at Bob (point C in Fig.l) within one gating 
window is 

Wpd) = N M o d VcHVDMu(N A %} + (N^ k )At + N Mod (Nj$ SRS ) (8) 

Here (N A A g E ), (N^J k } and (N^ SRS ) are determined from equations (3), (4) and (6) 
respectively. At the RHS of (8), the superscripts are used to refer to the location where 
noise photons are evaluated. 

The gain of EDFA is adjusted with the channel transmittance r] c h to maintain a 
constant P out (point B in Fig.l) of the classical channel. In this paper, we assume 
G = 100/77 cft . 

Using (3-8), we calculated the number of noise photons (Ng°p D ) as a function of the 
fiber length. We have assumed that P out = OdBm and G = 100/^. Other simulation 
parameters are summarized in Table 2. Fig. 3 shows the simulation results: at short 
distances, the main contribution of noise photons is the leakage from the classical 
channel; while at long distances, most of the noise photons are from SASRS. Noise 
contributed by the EDFA is negligible. 

In the BB84 QKD system, noise photons trigger random detection events and 
contribute to quantum bit errors. This contribution can be included in the total 
background count rate: 

Y = Y ° + VBob (N™ D ) (9) 
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Figure 3. The number of noise photons arrived within Ins gating window of SPD. 
The power of the classical signal output from the communication fiber is assumed to 
be P out = OdBm regardless of the fiber length. The EDFA gain G = 100 /r) c h- Other 
simulation parameters are summarized in Table 2. 



Here Yq is the background rate of the original QKD system. The error rate of 
background counts is assumed to be eo = 0.5. 

We assume that Alice and Bob perform perfect decoy state BB84 protocol with 
infinite decoy states. In the asymptotic limit of infinitely many signals sent by Alice, 
the secure key rate (per signal sent by Alice) is given by [H] 

R = \[Qx ~ f{E l x)QnH 2 {E il ) - Q x H 2 {e x )). (10) 

Here \i is the expected photon number of the signal state, Q^, are the gain and the 
overall quantum bit error rate (QBER) of signal states, while Qi, e\ are the gain and 
the QBER of single-photon components, / > 1 is the inefficiency factor of the error 
correction algorithm. The estimated values of the above parameters are given by [15] 



Q^ = Y + l-e-^ (11) 

E» = [e Y + e det (l - e'^/Q^ (12) 

Q 1 =(Y + Ti)ne-» (13) 

e a = (e Y + e det rj)fMe- fi /Q 1 (14) 



Here e^ is the probability that a single photon hits the wrong detector when Alice and 
Bob choose the same basis. r\ = r/chrjDMuVBob is the overall efficiency, where r\B b is the 
efficiency of Bob's system. 

Given the noise photons shown in Fig. 3, we calculated the secure key rate of the 
decoy BB84 QKD system using (9-14). The simulation parameters are chosen to be [TT] 
&det = 0.003, Yq = 5 x 10~ 6 , r] Bob = 0.038 and / = 1.22. The simulation results show 
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that no secure key can be generated at any distance. We remark that the simulation 
results are not sensitive to the actual values of e^t and F ° since the quantum bit errors 
are mainly contributed by the noise photons due to multiplexing. 

4-2. A homodyne detection based scheme: GMCS QKD 

The GMCS QKD has drawn a lot of attention for its potential high secure key rate, 
especially at relatively short distances [HI EHJ EU [32] . In this protocol, instead of 
performing single photon detection, Bob measures either the phase quadrature or the 
amplitude quadrature of a weak coherent state by using a homodyne detector. The 
strong LO used in homodyne detection acts as a "mode" filter: only noise photons in 
the same spatiotemporal and polarization mode as the LO can be detected, while noise 
photons in unmatched modes will be effectively suppressed. However, if the number 
of noise photons in unmatched modes is comparable to the photon number of the LO, 
their contributions cannot be neglected [T2] . 

Refer to Fig.l, in the GMCS QKD protocol, the optical power of each quantum 
channel is mainly determined by the operation rate and the average photon number of 
LO. Given a 1MHz operation rate and a LO of 10 8 photons per pulse, the optical power of 
each quantum channel is around — 19dBm, which is significantly lower than the average 
power of the classical signal (OdBm). Here, we simply neglect the crosstalk between 
quantum channels. If a pair of QKD users can access multiple quantum channels, the 
achievable secure key rate will scale with the channel number. We further assume that 
all the classical channels are at longer wavelengths of the QKD signals and the total noise 
contributed by all the classical channels is simply a summation of the noise contributed 
by each channel. 

4-2.1. Noise photons in matched mode Noise photons in the same spatiotemporal and 
polarization mode as the LO are contributed by in-band ASE and SASRS. The number 
of noise photons in matched mode arrived at Bob (point C in Fig.l) is 



Here, the factor 1/2 is due to the polarization selection of the LO, m is the number of 
classical channel, (N$ E ) and (Ng^ SRS ) are determined by (3) and (6). 

Both the ASE and the SASRS can be modeled as output from a chaotic source with 
Bose— Einstein photon statistics [331 US]- For a chaotic source with an average photon 
number of (n), the quadrature variances in shot noise unit are given by [35J: 



Using (15) and (16), the "excess noise" contributed by noise photons in matched 
mode is given by 



(Nqmcs) = lrn(VckVDMu(N { A %) + (N^srs)) 



(15) 



(AX) 2 = (AY) 2 = 2(n) + 1 



(16) 



Sin — 2VBob{NQ MCS ) 



(17) 
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4-2.2. Noise photons in unmatched modes In the normal working condition, the pulse 
width of the LO is significantly smaller than the integration time of the homodyne 
detector. The latter is determined by the bandwith Af of the homodyne detector and 
can be estimated by AT = 2 nAf Similar to (8), which gives the noise photon number 
within one gating window At of the SPD, the number of noise photons in "unmatched 
modes" measured at Bob (point C in Fig.l) within the time window of AT is given by 

AT 

(N° G tcs) = ^( N spd) (18) 

Again, we model the noise photons in unmatched modes as the output of a chaotic 
source. The photon number variance of a single mode chaotic light is given by [35J: 

(An) 2 = (n) 2 + (n) (19) 

where (n) is the average photon number. 

However, the integration time AT of the homodyne detector is normally much larger 
than the coherence time r c of the noise photon [36]. Under this condition, the photon 
number statistics follows Poisson distribution (37j. Thus the "excess noise" contributed 
by noise photons in unmatched modes is 

Eout = riBo b (N° G t C s)/(n LO ) (20) 

Note in (20), we have added in a factor of l/(n,Lo) to express e out in shot noise unit, 
where (tilo) is the average photon number of the LO. 

For example, if we assume that the bandwidth of the homodyne detector is 
Af = 1MHz, then AT = 0.16/xs. Since we have assumed that the gating window 
of the SPD is At = Ins, from (18), (Ng^ cs ) = 160 x (N%$ D ). Based on the results in 
Fig. 3, we expect that the total number of noise photons in unmatched modes (^gmcs) 
is in the order of 10 2 . If we further assume that (tilo) = 10 8 , then e out is in the order of 
10 -6 , which is negligible. In this paper, we will neglect the contribution from photons 
in unmatched modes. 



4-2.3. Secure key rate of the GMCS QKD Under the "realistic model" [17J, the secure 
key rate (per signal sent by Alice) of the GMCS QKD with "reverse reconciliation" 
protocol is given by [32] 

AI = 1 I AB - X be (21) 

where 7 is the efficiency of the reverse reconciliation algorithm, and [jj 

Iab = \ log 2 [(V + + Xtot)] (22) 

XBE = e(^ ) + e(^)-e(^)-e(^) (23) 

with Q(x) = (x + l)log 2 (x + 1) - xlog 2 x; Xtot = Xiine + Xhom/Vch] Xime = 1/Vch - 1 + e; 
Xhom = (1 + v el )/i - 1; a 2 2 = \{A ± ^ 2 - 4S); = \{C ± y/CP - AD); A = 

§ To avoid confusion, some symbols are different from the ones used in 325 
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2. fi V\ r B+Vch(V+Xiine)+Axho m . 
' ° Vch(V+Xtot) ' 

variance of the coherent state 



prepared by Alice, rf = r] DMU rj Bob is the equivalent efficiency of Bob's system, e denotes 
noise contribution from outside of Bob's system, which can be further separated into 
two terms: 



Here e is contributed by the original GMCS QKD system. e in is the excess noise due 
to multiplexing with the classical channels and can be determined from (17). Note that 
e in is determined at Bob's side, while e and e are referred to input. 

We have performed numerical simulations using parameters in Table 2 and the 
results are shown in Fig.4. In Fig.4, the secure key rates have been calculated under 
5 different conditions: (1) No classical signal; (2) One non-adjacent classical channel 
(OdBm); (3) One adjacent classical channel (OdBm); (4) 38 non-adjacent classical 
channels (OdBm per channel); (5) One non-adjacent classical channel (OdBm) with our 
100MHz homodyne detector (see details in Section 4.2.4). In all the above cases, the 
optical power is defined at the output of the communication channel (point B in Fig.l). 
Note under condition (3), we have assumed that the isolation of MUX/DEMUX between 
adjacent channels is — 40dB, which is achievable with commercial products [JJ. From 
Fig.4, it is possible to multiplex the GMCS QKD with a OdBm classical channel without 
significantly reducing its performance. Even multiplexed with 38 classical channels 
(OdBm power each channel), the GMCS QKD still has a secure distance around 10km. 
We remark that Eq. (23) was derived from the realistic model [17J, where Eve cannot 
take advantages of noise contributed by Bob's system [38]. In the more conservative 
"general model" [17], where Eve can control losses and noise in Bob's system, the secure 
distance is much shorter [31]. In practice, the realistic model has been commonly used 



4-2.4- Preliminary experimental results with a 100MHz homodyne detector Recently, 
we have developed a 100MHz bandwidth shot-noise limited optical homodyne detector 
[39] . With a LO photon number of 4 x 10 8 , the electrical noise of the homodyne detector 
is IOgLB below the shot noise. 

The noise photons output from the DWDM system (point C in Fig.l) has been 
fed into the 100MHz homodyne detector. The fiber length used in this experiment is 
20km. The variance of the output of the homodyne detector has been measured as 
a function of the LO power under three conditions: (1) No classical channel; (2) One 
classical channel (P ou t — 3AdBm, G = lOO/r^) is placed at channel 2, channel 8 is used 
as the quantum channel; (3) One classical channel (Pout = 3AdBm,G = 100/r] c h) is 
placed at channel 2, channel 1 is used as the quantum channel. Note under condition 



e = e + 



(24) 



VchVDMUVBob 



[261 EH EE]. 



|| For the specific device we have tested, the isolation between adjacent channels varies from channel 
to channel in the range of — 30dB to —50dB 
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Fiber length (km) 



Figure 4. The secure key rate of the GMCS QKD under the "realistic model" . The 
secure key rates have been calculated under 5 different conditions: (1) No classical 
signal; (2) One non-adjacent classical channel (OdBm); (3) One adjacent classical 
channel (OdBm); (4) 38 non-adjacent classical channels (OdBm per channel); (5) 
One non-adjacent classical channel (OdBm) with our 100MHz homodyne detector 
(see details in Section 4.2.4). Note under condition (3), we have assumed that the 
isolation of MUX/DEMUX between adjacent channels is — AOdB. Other parameters 
are summarized in Table II. 



(3), the quantum channel is placed at the adjacent channel of the classical channel. 
Under both condition (2) and (3), the additional excess noise due to multiplexing has 
been estimated from (17) to be less than 0.01 (in shot noise unit). On the other hand, 
the measurement uncertainty of our homodyne detection system on determining the 
noise variance has been measured to be <j = 0.024 (in shot noise unit) 0. As shown in 
Fig. 5, there is no observable difference among the 3 measurement results. To determine 
the magnitude of multiplexing noise more accurately, the measurement accuracy of the 
homodyne detector has to be further improved. 

This raises up a question about how to apply the realistic model using this 
homodyne detector, since Bob has to estimate the excess noise from his system and the 
one from outside separately. To determine the excess noise contribution from outside 
of Bob's system (e), one common way is to subtract the total observed noise by the 
vacuum noise associated with the transmission efficiency and the electrical noise (v e i) 
contributed by Bob's homodyne detector. Obviously, the minimal e resolvable by this 
method is determined by the measurement uncertainty <; of the homodyne detector El- 
Note that the simulation results shown in Fig. 4 (except the one acquired under 
condition (5)) are based on parameters of the GMCS QKD system from [32], where a 
low speed (1MHz bandwidth) and low electronic noise (20dB below the shot noise) 



% The measurement uncertainty is defined as three times of the standard deviation 

+ Assume that the transmission efficiency (thus the vacuum noise) can be determined accurately 
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-* — No classical signal 

- — Ch8/Ch2: Quantum/Classical 

-e — Ch1/Ch2: Quantum/Classical 




Power of LO (mw) 



Figure 5. Experimental results with a 100MHz homodyne detector. The total 
variance of the output of the homodyne detector has been measured as a function 
of the LO power under three conditions: (1) No classical channel; (2) One classical 
channel (P ou t = 3AdBm, G = 100/r) c h) is placed at channel 2, channel 8 is used as the 
quantum channel; (3) One classical channel (P ou t — 3AdBm,G — 100/r) c h) is placed 
at channel 2, channel 1 is used as the quantum channel. The fiber length is 20km. 
From the experimental results, there is no observable difference among the 3 curves. 

homodyne detector was employed. So we have assumed that the measurement 
uncertainty of the homodyne detector is small enough to resolve the excess noise due 
to multiplexing. In our preliminary experiment, however, the homodyne detector has 
a larger bandwidth (100MHz) but also higher electronic noise (lOdB below the shot 
noise). As a conservative estimation of e, we could replace (24) by 

6 = e + ^ + - (25) 

VchVDMUVBob VchVDMUVBob 

Using (25), we have performed numerical simulations based on parameters of the 
100MHz homodyne detector: v e \ = 0.1 and q = 0.024. Other parameters are shown 
in Table 2. We have assumed that one non-adjacent classical channel (OdBm) is 
multiplexed with the quantum channel. The simulation result is also shown in Fig. 4 
(under condition (5)) where the secure distance is about 14km. Note that the secure 
key rate under condition (5) is significantly lower than that under condition (2). This is 
mainly due to the larger measurement uncertainty of the 100MHz homodyne detector. 
To apply the realistic model more efficiently, we may need a more accurate way than 
the one given by (25) to estimate e. 

We would like to end this section with a few comments on the realistic model 
adopted in the GMCS QKD. In all the security proofs mentioned in Section 1, one 
underlying assumption is that both Alice and Bob's QKD systems are fabricated by 
trusted vendors and these devices are placed inside Alice and Bob's local secure stations 
which cannot be accessed by Eve. So it might be reasonable to assume that Eve cannot 
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control the internal parameters of Bob's system (the realistic model). However, to 
justify the above assumption in practice, we may need to develop special techniques to 
estimate each system parameter accurately without compromising the security of the 
QKD system. As we have discussed in Section 1, to apply the security proof of an 
idealized QKD protocol to a practical QKD system, all the underlying assumptions and 
implementation details have to be studied carefully. 

5. Conclusion 

In summary, we have studied the feasibility of conducting QKD together with classical 
communication through the same fiber by employing C-band DWDM technology. The 
impact of the classical channel to the quantum channel has been investigated for both 
QKD based on single photon detection and QKD based on homodyne detection. Our 
studies show that the latter can tolerate a much higher level of contamination from the 
classical channel than the former. We have performed simulations based on both the 
decoy BB84 QKD protocol and the GMCS QKD protocol. With commercial DWDM 
components, our simulation results show that it is possible to multiplex the GMCS 
QKD with a OdBm classical channel without significantly reducing its performance. 
Even multiplexed with 38 classical channels (OdBm power each channel), the GMCS 
QKD still has a secure distance around 10km. 

Although the LO is assumed to be a single mode coherent state in this paper (which 
is not difficult to achieve in practice), it doesn't have to be transform-limited. 

The noise photons in the BB84 QKD system could be further reduced by employing 
narrow spectral filter and temporal filter. For example, in a 45pm spectral filter 
has been employed to further cut off noise. In [40], the author suggested to use an 
additional optical gate to further suppress noise photons in time domain. In principle, 
it is possible to selectively detect photons in only one spatiotemporal mode by using an 
optimal combination of spectral and temporal filters [41]. However, in practice, both the 
spectral filter with extremely narrow bandwidth and the temporal filter with extremely 
narrow time window are difficult to fabricate and lossy as well as unstable (subject 
to minute changes in temperature, pressure, etc.). Furthermore, Alice's signal has to 
be transform-limited to pass through these filters effectively. This also requires careful 
dispersion management in the communication channel. 

Financial support from CFI, CIPI, the CRC program, CIFAR, MITACS, NSERC, 
OIT, and Quantum Works is gratefully acknowledged. 
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